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\ Abstract 

o . 

. In it was proposed, that the nearest-neighbor distribution P{s) of the spectrum of the 

' Bohr-Mottelson model is similar to the semi-Poisson distribution. We show however, that 

. P{s) of this model differs considerably in many aspects from semi-Poisson. In addition we 

give an asymptotic formula for P(s) as s — > 0, which gives P'(0) = 7rV3/2 for the slope 
^ . at s = 0. This is different not only from the GOE case but also from the semi-Poisson 

prediction that leads to /"(O) = 4 



Q , The motivation for this comment stems from the impression, that the article "models of 

intermediate spectral statistics" can easily be misinterpreted in two ways: One might be led 
to believe that (i) the semi-Poisson distribution is universal, and (ii) the universality class of 
"intermediate statistics" is as well defined and established as for example the Poisson ensemble 
\ or the GOE. In this comment, we will argue, that both statements are wrong. 

The purpose of [P is to present models which could constitute a "third" universality class of 
systems which show so called "intermediate statistics", previously introduced by Shklovskii in 
. The Poissonian and the Gaussian ensemble (for definiteness, consider orthogonal ensembles 
only) are considered as the first two universality classes in this list. 

As in the Poissonian and in the GOE case, where the respective members have common and 
unique statistical properties, one would expect the same to hold for the models with intermediate 
statistics. In Q the authors concentrate on the distribution of nearest neighbor spacings. In 
the Poissonian case it is given by P{s) = exp(— s), in the GOE case it is close to the well known 
Wigner surmise P{s) ~ (7r/2) exp(— 7r/4s^), whereas in the case of the "intermediate statistics" 
the candidate proposed in [|| is the semi-Poisson distribution P{s) = 4s exp(— 2s). 

In what follows we will show, that the level spacing distribution in the case of the Bohr- 
Mottelson model, that represents one of the candidate systems for the intermediate statistics 
mentioned in Q is in fact very different from the proposed semi-Poisson distribution. This 
discrepancy can be found in fact already on a figure published in Q but without discussing the 
problem. However, |Q] gives an overview over the statistical properties of different variants of 
the Bohr-Mottelson model. 
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Figure 1: (a) P{s) for the model (|l|) compared with the the semi-Poisson distribution, (b) the 
same as in (a) for short distances. In addition the theoretical result (^) is drawn as a dashed 
line and the GOE result as a dashed dotted curve. 



In the following matrix model (originally introduced by Bohr and Mottelson) is presented 
as a possible candidate showing statistical properties similar to the semi-Poisson distribution: 

Hmn — ^n^mn ~\~ tmtn • (1) 

H is a N X A^-matrix, e„ are mutually independent random variables uniformly distributed over 
a finite interval, and the t„ are chosen with equal absolute value squared = r. The authors 
of HI sketch a procedure for calculating analytically the 2-point correlation function. For small 
distances it should agree with P{s), so that one can derive the slope P'(0) = 7r\/3/2 of the 
spacing distribution at s = 0. This slope is different from the GOE case, where it equals to 7r^/6 
as well as from the slope of the semi Poisson distribution that equals to 4. The difference to the 
GOE case is remarked, but a similar remark on the difference to the semi-Poisson is avoided. 
An unprejudiced reader might believe, that the correlation properties of the matrix model are 
similar to semi-Poisson, even though that the spectral statistics of this model differ remarkably. 

In Figure 1 we present the numerical result for P{s) obtained for an ensemble of 1000 matrices 
of dimension 750. For the statistical analysis we used only one third of the states in the center 
of the spectral region. The numbers are uniformly distributed over an interval [—1,1] and the 
elements of the vector t are chosen as ti = y/a/JjTpN), where p is the level density in the center 
of the spectrum, A'^ is the dimension of the matrix and a = 10 is the coupling constant. (We 
checked, that larger coupling does not change the numerical results). Figure l.a demonstrates 
the qualitative differences in the behavior of P{s) between the random matrix model and the 
semi-Poisson distribution. The slope at s = is smaller, the maximum of P{s) is slightly shifted 
to the right, and for values s > 1 it shows significant deviations well above of the statistical 
error. Hence the spacing distribution of the Bohr- Mottelson model is not close to semi-Poison. 

Figure l.b. shows a magnification of the interval < s < 1/2 using the same data as in 
figure l.a. In this figure we additionally plotted the asymptotic result (^ for the present model 
as a dashed line. The basic idea of the approximation is the following: In order to get a short 
distance between two neighbored levels in the spectrum of three eigenvalues of Hq have to 
come close together. Then the levels which are farther away can be neglected. Therefore we can 
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restrict the sum 

K{E) = y: 



1=1 



whose roots define the eigenvalues of H, to those terms with the three consecutive eigenvalues. 
Resolving for the two roots, calculating their distance, and averaging over the levels ej leads to 
the following formula 



Tv/2 

9s f , exp 



3s 



(cos (j) + sin 0)/v/l + sin(2(/))/2 



Pis) = -r / # — ^-7-7TT7wT^ ■ (3) 



4: J ^ l + sin(2</>)/2 



The dashed curve in Figure l.b. is obtained from a numerical integration of At short 
distances, this approximation describes the numerical data much better than the semi-Poisson. 
A Taylor expansion of the integrand of @ gives -P'(O) = 7r^/3/2 for the slope at s = 0, the same 
result as found in |||]. 

A detailed numerical investigation of several statistical properties of the type of models can 
be found also in |^, ^, 

To conclude, if one wants to insist on the introduction of a "third" universal ensemble, one 
should possibly use a criterium similar to the following (cited from [Q: "... the main features are 
(i) the existence of level repulsion (as in random matrix ensembles), and (ii) slow (approximately 
exponential) fall-off . . . " . Even though this definition is quite "spongy" , it seems to be the only 
way to make sure, that the systems discussed fit into this class. 
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